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ABSTRACT 

We present strong detections of methyl cyanide (CH 3 CN), vinyl cyanide 
(CH 2 CHCN), ethyl cyanide (CH 3 CH 2 CN) and cyanodiacetylene (HC 4 CN) 
molecules with the Green Bank Telescope (GBT) toward the Sgr B2(N) molec- 
ular cloud. Attempts to detect the corresponding isocyanide isomers were only 
successful in the case of methyl isocyanide (CH 3 NC) for its Jk = 1q — Oo tran- 
sition, which is the first interstellar report of this line. To determine the spatial 
distribution of CH 3 NC, we used archival Berkeley-Illinois-Maryland Association 
(BIMA) array data for the Jk = ~ 3 k {K = — 3) transitions but no emission 
was detected. From ab initio calculations, the bonding energy difference between 
the cyanide and isocyanide molecules is >8500 cm -1 (>12,000 K). Thus, cyanides 
are the more stable isomers and would likely be formed more preferentially over 
their isocyanide counterparts. That we detect CH 3 NC emission with a single 
antenna (Gaussian beamsize(f2s)=1723 arcsec 2 ) but not with an interferometer 
(f2s=192 arcsec 2 ), strongly suggests that CH 3 NC has a widespread spatial distri- 
bution toward the Sgr B2(N) region. Other investigators have shown that CH 3 CN 
is present both in the LMH hot core of Sgr B2(N) and in the surrounding medium, 
while we have shown that CH 3 NC appears to be deficient in the LMH hot core. 
Thus, large-scale, non-thermal processes in the surrounding medium may account 
for the conversion of CH 3 CN to CH 3 NC while the LMH hot core, which is dom- 
inated by thermal processes, does not produce a significant amount of CH 3 NC. 
Ice analog experiments by other investigators have shown that radiation bom- 
bardment of CH 3 CN can produce CH 3 NC, thus supporting our observations. We 
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conclude that isomers separated by such large bonding energy differences are dis- 
tributed in different interstellar environments, making the evaluation of column 
density ratios between such isomers irrelevant unless it can be independently 
shown that these species are co-spatial. 

Subject headings: ISM: abundances - ISM: clouds - ISM: individual (Sagittarius 
B2(N)) - ISM: molecules - radio lines: ISM 



1. INTRODUCTION 

The observations of isomers in the interstellar medium is becoming increasingly more 
common as new and larger molecules are being detected. Currently, 135 molecules have been 
detected in interstellar and circumstellar environments. Excluding diatomic molecules, ~30% 
of all interstellar molecules have observed isomeric counterparts. In addition, as the number 
of atoms in an interstellar molecule increases, so do the number of isomers. For example, 
of the 9 interstellar molecules containing 8 atoms, 5 are isomers. With an empirical for- 
mula C2H4O2, acetic acid (CH3COOH), glycolaldehyde (CH 2 OHCHO) and methyl formate 
(CH3OCHO) are isomers; and for H 2 C 6 , hexapentaenylidene and triacetylene (HC 6 H) are 
isomers. Isomerism may provide important clues to the formation of interstellar molecules. 
As such, an important class of isomers are cyanide (-CN) and isocyanide (-NC) species. 

Hydrogen cyanide (HCN), methyl cyanide (CH 3 CN), vinyl cyanide (CH 2 CHCN), and 
ethyl cyanide (CH 3 CH 2 CN) all have large abundances and are easily detectable in a num- 
ber of interstellar environments and are readily seen in the spectra of several molecular line 
surveys (e.g. Friedel et al. 2004, Nummelin et al. 1998). However, while there have been 
detections of HNC in interstellar and cometary environments, there have been no confirmed 
detection of the isocyanides CH 3 NC, CH 2 CHNC, and CH 3 CH 2 NC The first tentative detec- 
tion of CH 3 NC was toward Sgr B2(OH) by Cernicharo et al. (1988) from observations of the 
J = 4 — 3, J = 5 — 4, and J = 7 — 6 transitions. Though the lines were highly contaminated 
from other known interstellar species, they calculated an CH 3 NC/CH 3 CN abundance ratio 
of ~0. 03-0. 05. Since this initial claim, there have been no further observations to confirm 
the existence of interstellar CH 3 NC, in particular toward the LMH hot core of Sgr B2(N), 
the best known source of large molecules. Thus, in this work, we confirmed the existence 
of CH 3 NC and attempted to observe other large isocyanides by observing the same (or sim- 
ilar) low energy transitions of the cyanides CH 3 CN, CH 2 CHCN, CH 3 CH 2 CN, and HC 4 CN 
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and their corresponding isocyanide isomers toward the the high-mass star forming region, 
Sgr B2(N). 

2. OBSERVATIONS AND RESULTS 

The observations using the NRAO 100-m Robert C. Byrd Green Bank Telescope (GBT) 
were made in 2004 February 24, February 29 and March 29. Table 1 lists the molecular line 
parameters of the transitions of the cyanides and isocyanides sought. The species, calculated 
rest frequencies, transition quantum numbers, line strengths, upper level energies, and con- 
tinuum levels are listed in the first six columns of Table 1. The GBT half-power beamwidth 
(8b) is given as 740"/ v [GHz]. The K-band receiver set used has a frequency range covering 18 
GHz to 22.5 GHz. For the observations reported herein, the average beamwidth < 9b >=39" 
corresponds to an average Gaussian beamsize of < Qb > =71 "^b/[4 ln(2)] = 1723 arcsec 2 and 
an average beam efficiency < r] B >~0.8. The GBT spectrometer was configured in its 8 
intermediate frequency (IF), 200 MHz, 3- level mode which permits observing four 200 MHz 
frequency bands at a time in two polarizations through the use of offset oscillators in the IF. 
This mode affords 24.4 kHz channel separation. Antenna temperatures are on the T* A scale 
(Ulich & Haas 1976) with estimated 20% uncertainties. The Sgr B2(N-LMH) J2000 pointing 
position employed was a— 17 ft 47 m 19 s .8, 5 = -28°22'17" and the LMH systemic LSR source 
velocity of +64 km s _1 was assumed. Automatically updated dynamic pointing and focusing 
corrections were employed based on realtime temperature measurements of the structure 
input to a thermal model of the GBT; zero points were adjusted typically every two hours 
or less using the calibrators 1626-298 and/or 1733-130. Data were taken in the OFF-ON 
position-switching mode, with the OFF position 60' east in azimuth with respect to the ON 
source position. A single observation consisted of two minutes in the OFF source position 
followed by two minutes in the ON source position. The two polarization outputs from the 
spectrometer were averaged in the final data reduction process to improve the signal-to-noise 
ratio. 

Figure 1 shows our attempt to detect cyanide and isocyanide species toward Sgr B2(N) 
with the GBT. The panels in the left column show the detections of CH 3 CN, CH 2 CHCN, 
CH3CH2CN and HC4CN. We have successfully detected all species in the left column and for 
the first time the JR a ,K c — 2n — lio transition of interstellar CH 3 CH 2 CN is reported. In the 
panels containing the transitions of CH 3 CN, CH 2 CHCN and CH 3 CH 2 CN, the intensity scale 
(ordinate) decreases, showing that the emission from the larger species is in general weaker. 
The abscissa is the radial velocity with respect to the local standard of rest (LSR) calculated 
for the rest frequency of the transition shown at the top left of each panel assuming a source 
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velocity of +64 km s 1 . Dashed lines show LSR velocities of +64, +73 and +82 km s 1 . 
Molecules with no hyperfine (HF) structure previously observed with the GBT such as 
propenal (CH 2 CHCHO), propanal (CH 3 CH 2 CHO), and glycolaldehyde (CH 2 OHCHO) show 
absorption components at +64 and +82 km s _1 and an emission component at +73 km s _1 
(Hollis et al. 2004a, b). In general, the cyanide spectra shown in the left panel of Figure 1 
are consistent with this source velocity structure. In the cases of CH 3 CN and CH 2 CHCN, 
however, the source velocity structure is confused by molecular HF splitting. The fiducial 
on each HF component denotes the position corresponding to +64 km s -1 . Thus, it can be 
seen that only the F=0-1 components are clear of any complication. The panels in the right 
column of Figure 1 show the passbands containing the transitions of CH 3 NC, CH 2 CHNC, 
CH 3 CH 2 NC and HC 4 NC. We clearly detect the absorption and emission components at +64, 
+73 and +82 km s _1 from CH 3 NC and there is possible narrow emission from CH 2 CHNC 
near +64 km s _1 . The primary absorption component of CH 3 NC at +64 km s _1 is more 
than an order of magnitude weaker in intensity than the corresponding feature of CH 3 CN, 
however this is the first interstellar report of this line. No emission or absorption is seen 
from either CH 3 CH 2 NC or HC 4 NC beyond a 1 a detection limit (~ 2 mK). 

The Jk = 4:r~^k {K = 0—3) CH 3 NC observations using the Berkeley-Illinois-Maryland 
Association (BIMA) Millimeter Array 1 were made in 2001 April 20 in its C configuration to- 
ward Sgr B2(N-LMH). The pointing position was the same as used for the GBT. Table 2 lists 
the molecular line parameters of each species in the CH 3 NC passband. The species, the calcu- 
lated rest frequencies, the transition quantum numbers, the product of the line strength and 
the square of the dipole moment (S/i 2 ), and the upper level energies are listed in the first five 
columns of Table 2. The synthesized beamwidth for these observations is 8 a x0b = 23."9x7."l 
corresponding to a Gaussian beamsize of f2#=192 arcsec 2 . The quasar 1733-130 was used to 
calibrate the antenna based gains and Uranus was used as the flux density calibrator. The ab- 
solute amplitude calibration of these sources is accurate to within ~20%. The passbands were 
calibrated automatically during data acquisition; a technical description of auto-calibration 
can be found at http://astron.berkeley.edu/~plambeck/technical.html. The spectral win- 
dow containing the CH 3 NC transitions has a bandwidth of 50 MHz and is divided into 128 
channels for a spectral resolution of 0.39 MHz per channel. All data were combined, imaged, 
and self-calibrated using the MIRIAD software package (Sault, Teuben & Wright 1995). 

Figure 2 shows fiducials for the Jk = — 3k (K — — 3) transitions of CH 3 NC in a 
passband dominated by strong emission from CH 3 CH 2 CN which is largely confined to the 
LMH hot core (Liu & Snyder 1999). The rest frequency located in the upper left of the panel 
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corresponds to the Jx a ,K c = 928 — 827 transition of CH 3 CH 2 CN for an assumed LSR velocity 
of +64 km s -1 . Dashed lines denote LSR velocities of +64, +73 and +82 km s -1 . The 
transition of CH 3 CH 2 CN clearly shows strong emission components at the LMH systemic 
velocity of +64 km s -1 and also at +73 km s _1 . The J K = A K - 3 K (K = - 3) CH 3 NC 
transitions are labeled in this passband assuming an LSR velocity of +64 km s" 1 and, except 
for the Jk = 4 3 — 3 3 transition, all lie in a region that is clear from any other molecular line 
emission. However, no transition of CH 3 NC is detected above a 1 a detection limit (~0.3 Jy 
beam -1 ). The 1 a rms noise level fiducial is shown at the left of the panel. 



We have presented two disparate sets of observations: GBT single antenna results with 
a Gaussian beamsize of 0^=1723 arcsec 2 , and BIMA array results obtained with a Gaussian 
beamsize of f2 B =192 arcsec 2 . Given the GBT beamsize is a factor of 9 greater than the 
BIMA beamsize, the spatial scale of the GBT will couple better to the extended emission 
toward Sgr B2(N). Conversely, the BIMA array will tend to couple better to the 28 arcsec 2 
(~5" diameter) of the LMH hot core. Hence, we expect that the GBT observations will 
be characterized by a lower excitation temperature (T ex ) than the excitation temperature 
associated with the BIMA observations. First, we will derive the cyanide and isocyanide 
column densities from the GBT observations. Then, we will show from BIMA observations 
that CH 3 NC must be extended. 

For the GBT observations, following the formalism outlined in Hollis et al. (2004b), the 
total beam-averaged column density of a molecular species obtained with a single element 
radio telescope is given as: 

For emission: 



3. 
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In both equations, the line shapes are assumed to be Gaussian; rjB is the telescope 
beam efficiency; T ex is the excitation temperature; AT\AV is the product of the fitted line 
intensity (mK) and line width (km s -1 ); Q r is the rotational partition function; Sfi 2 is the 
product of the transition line strength and the square of the dipole moment (Debye 2 ); and 
E u is the upper rotational energy level (K). In equation (1), v is the transition frequency 
(MHz) and Tb g ~2.7 K is the cosmic background temperature. In equation (2), T c is the 
source continuum temperature (K), and E x is the lower rotational energy level (K). 

Using the GBT at K-band, Hollis et al. (2004b) found a glycolaldehyde excitation tem- 
perature of T ex =8 K that satisfied the detected integrated line intensities of the absorption 
components at +64 and +82 km s _1 and the emission component at +73 km s _1 . Since the 
LSR velocity components of the cyanide and isocyanide molecules in Figure 1 appear to be 
similar, we used T ex =8 K to find the total column density of each of the detected species. 
Columns (7) to (12) of Table 1 summarize the Gaussian fitting intensities and line widths 
for each detected transition or 3 a upper limits for the intensities of non-detections at LSR 
velocities of +64, +73 and +82 km s _1 . 

In order to determine the intensity and line width of the individual HF and velocity 
components of the CH 3 CN and CH 2 CHCN lines, multiple Gaussian fits were used. By fixing 
the frequency separation between HF components, we fit the +64, +73 and +82 km s -1 
velocity components for each HF component. If the fitting routine found a solution for the 
intensity of a velocity component less than the 3 a rms noise level, it was omitted, and the fit 
recalculated for one less velocity component. The results are listed in columns 7-12 of Table 
1. In Figure la, the solid line is the result of Gaussian fitting the 3 HF structure components 
to the +64, +73 and +82 km s _1 LSR velocity structure (see Table 1). Furthermore, as a 
residual to this methodology, the small absorption dip near 50 km s -1 in Figure la is not 
due to a HF component nor is it a velocity component of CH 3 CN. 

The column densities listed in column (13) of Table 1 are the sum of all detected 
velocity components, or, in the cases of CH 3 CN and CH 2 CHCN, of all velocity and HF 
components. The column densities of the cyanide molecules are within a factor of 3 or 4 of the 
CH 3 CN column density of (~10 14 cm -2 ). Furthermore, the column density measured from 
the +64 km s -1 velocity component accounts for ~75% of the total CH 3 CN column density. 
The +64 km s _1 velocity component is often the only one resolved at higher frequencies 
and is the basis for the majority of CH 3 CN column density measurements at 1 and 3 mm 
wavelengths. 

The similarity in column densities between cyanides is also seen at higher resolution 
toward other sources (Remijan et al. 2004b). The column density ratio CH 3 NC/CH 3 CN is 
~0.02 and is relevant if and only if the two species are co-spatial. The possible detection of 
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CH 2 CHNC gives a column density ratio CH 2 CHNC/CH 2 CHCN of -0.005. 

From Remijan et al. (2003), the total beam averaged column density obtained with an 
interferometer is: 

fAIdvQ r e^ T ^ 2n 2 

Nt = 2Mx Be a e h ^ xl ° cm < 3 > 

In equation (3), J Aldv is the product of the fitted line intensity (Jy beam" 1 ) and line width 
(km s -1 ), 9 a and 9 b are the FWHM beam widths (arcsec) and B = Q s /[^b + &s] is the 
beam filling factor where f2 B is the solid angle subtended by the synthesized beam of the 
interferometer and ils is the solid angle subtended by the source emission. For the BIMA 
array data, 9 a x9 b = 23."9x7."l. Columns (7) to (12) of Table 2 summarize the Gaussian 
fitting intensities and line widths for CH 3 CH 2 CN and the 3 a upper limits for the intensities 
of the non-detections at LSR velocities of +64, +73 and +82 km s -1 . The column density 
of CH 3 CH 2 CN in column (13) is based on the excitation temperature of methanol (CH 3 OH) 
of 170 K (Pei, Liu & Snyder 2000). 

No emission for the J = 4 — 3 (K = — 3) transition of CH3NC was seen in the 
BIMA array SgrB2(N) data. However, if CH 3 NC and CH 3 CN are co-spatial (i.e., our GBT 
CH3NC/CH3CN ratio of ~0.02 is valid), then observations of CH 3 CN by other investigators 
which yield large CH 3 CN column densities predict that we should have detected CH 3 NC 
with the BIMA Array if the source is compact. In an extensive survey encompassing the 
J = 12 — > 11, J = 13 — > 12 and J = 14 — > 13 transitions of CH 3 CN at 1 mm wavelengths, 
Nummelin et al. (2000) obtained a CH 3 CN beam averaged column density of ~4.2x 10 15 cm -2 
for a 352 K LMH hot core temperature. These observations were performed using the 
15m SEST telescope with a ~600 arcsec 2 beamsize. Owing to the lower spatial resolution 
(~2 km s" 1 channel -1 ) at 1 mm wavelengths and the large line widths (15-20 km s _1 ) of 
CH 3 CN lines toward Sgr B2(N-LMH) (Nummelin et al. 2000), it is not possible to resolve 
the individual velocity components of the CH 3 CN lines. Thus, the calculated column density 
is based solely on the +64 km s _1 LSR velocity component. In addition to the the column 
density and temperature, Nummelin et al. also found a beam filling factor of the CH 3 CN 
emission toward the Sgr B2(N) hot core of 0.014. In their data analysis, the column density, 
rotational temperature and beam filling factor were given as free parameters. From these 
values, an integrated line intensity was calculated for each transition and was compared 
to the observed integrated line intensity until a "best-fit" model was obtained. The best- 
fit model was found by finding the minimum of the reduced x 2 function (Nummelin et al. 
2000). Assuming the model explains the physics accurately, the measured values of the 
rotational temperature, column density and beam filling factor by Nummelin at al. toward 
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the Sgr B2(N) hot core are reported at a 68.3% confidence interval (1 a). Given the ~600 
arcsec 2 beamsize of the 15m SEST telescope, and a beam filling factor of 0.014, the measured 
source size of the CH 3 CN emitting region toward the Sgr B2(N) hot core is ~8.5 arcsec 2 . 
This size is very similar to the CH 3 CN cores toward W51el (7.9(6)arcsec 2 ) and W51e2 
(6.0(9)arcsec 2 ) (Remijan et al. 2004a). Furthermore, the size is also similar to the emitting 
regions of CH 2 CHCN and CH 3 CH 2 CN towards the Sgr B2(N) hot core measured by Liu & 
Snyder (1999) at subarcsecond resolution. Thus, ~8.5 arcsec 2 is a reasonable measure of 
the CH 3 CN emitting region toward the Sgr B2(N) hot core. Using the measured rotational 
temperature of 352 K, the scaled Nummelin et al. (2000) CH 3 CN column density within the 
192 arcsec 2 beamsize of the BIMA array is ~9.9xl0 16 cm -2 . Using the column density ratio 
of 0.02, the expected CH 3 NC column density of the +64 km s -1 LSR velocity component 
should be ~2xl0 15 cm -2 for co-spatial isomers. 

Solving equation (3) for J Aldv assuming a column density of ~2x 10 15 cm -2 , we should 
have detected a CH 3 NC integrated line intensity of ~57 Jy beam -1 km s -1 for the Jk = 
4o — 3o transition. Assuming a typical line width of 15 km s -1 for the LMH hot core 
of Sgr B2(N), the expected peak intensity of the 4o — 3o transition is ~3.8 Jy beam -1 . 
However, the BIMA array data shows no CH 3 NC emission greater than the 1 a rms upper 
limit of ~0.3 Jy beam -1 . The most likely explanation for the lack of a CH 3 NC detection 
is that the emission is extended and the array cannot couple to it. In the case of extended 
emission, the ratio of the measured line intensity, or in our case the 1 a rms upper limit of 
~0.3 Jy beam -1 , to the expected line intensity of ~3.8 Jy beam -1 is equal to the ratio of 
the solid angle subtended by the synthesized beam of the interferometer (VL B ) to the solid 
angle subtended by the source emission (fls)- This is because the synthesized beam of the 
interferometer is only collecting a fraction of the flux from the extended emission region which 
we assume to be uniformly distributed. Therefore our observations indicate that Qb/Qs < 
0.3/3.8. Solving for the source size assuming a Gaussian beamsize of 192 arcsec 2 , we find 
fis>2432 arcsec 2 or equivalently, a linear size of >100"(1.7'), which is much larger than the 
typical LMH hot core diameter of ~5" (Hollis et al. 2003). This strongly suggests that the 
spatial distribution of CH 3 NC is extended and would not show any significant concentration 
toward the LMH hot core. 

DeFrees, McLean, & Herbst (1985) theoretically calculated a CH 3 NC/CH 3 CN abun- 
dance ratio for dense interstellar clouds in the range of 0.1-0.4 by assuming that both species 
are formed from ion precursors and destroyed by gas phase, ion molecule reactions. The 
DeFrees et al. calculations show that when CH 3 + collides with HCN, a collision complex is 
formed that rapidly equilibrates to protonated methyl isocyanide (CH 3 NCH + ) and proto- 
nated methyl cyanide (CH 3 CNH + ) due to the large amount of energy produced from the 
initial collision. As the ensemble of ions continues to relax, the internal energy of CH 3 CNH + 
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becomes insufficient to overcome the barrier to isomerize to CH 3 NCH + . The end result finds 
CH 3 CN formed preferentially over CH 3 NC We performed ab initio calculations in order to 
quantify the energy differences between cyanide and isocyanide isomers in this work. Table 
3 lists each molecular species (col. 1), the zero-point corrected bonding energy (ZPCBE in 
col. 2) and the relative ZPCBE (col. 3). Column (2) ZPCBE represents the complete disso- 
ciation of the molecule into its constituent atoms. Note that the relative bonding energies 
between isomers are very large (>8500 cm -1 ). Thus, the cyanide molecules are much more 
stable than their isocyanide isomers, and therefore it is highly unlikely that any thermal 
process can satisfactorily explain the formation of isocyanide molecules. This indicates that 
a non-thermal process is necessary to convert a cyanide into an isocyanide species. 

In recent laboratory experiments performed by Hudson & Moore (2004), cyanide species 
were bombarded by both protons and UV-photons to determine what would happen to the 
structure of the molecule. Bombarding a pure ice sample of CH 3 CN at 24 K, the proton 
experiment yielded ketenimine (CH2CNH), another isomer of CH3CN, as well as CH3NC in 
a slightly lower abundance. Similarly, the UV-photon experiment at 12 K yielded CH3NC 
as the primary product. These products were also seen when CH3CN was embedded in N2 
ice and then irradiated. Hudson & Moore (2004) also note that the temperature in both 
experiments did not affect the formation of the end products. In an UV-photon experiment 
performed at 100 K, the destruction of CH 3 CN was faster than at 12 K, but the same end 
products were seen. This suggests that the reactions were driven more by the radiation 
dose than the sample temperature (Hudson & Moore 2004), reinforcing our suggestion that 
non-thermal processes may be the primary route to the formation of interstellar isocyanides. 
Laboratory experiments involving CH2CHCN, CH3CH2CN, and cyanoacetylene (HC2CN), 
all showed similar results for isocyanide isomers. It is likely that a similar reaction involving 
HC4CN will produce HC4NC. Furthermore, these laboratory experiments also suggest that 
all isocyanides are likely to be formed under the conditions which interstellar CH 3 NC was 
detected. Based on the column density ratios and upper limits in Table 1, we did not de- 
tect any emission from CH 3 CH 2 NC and HC 4 NC because of the possibility that we did not 
integrate long enough to detect these species. On the other hand, Hudson & Moore (2004) 
also conducted experiments with cyanides embedded in H 2 ice and found that hydrogen 
addition and subtraction reactions were taking place in both radiation bombardment reac- 
tions in which no isocyanides were formed. In these experiments HC2CN hydrogenated to 
CH2CHCN and then to CH 3 CH2CN. Based on the aforementioned laboratory results, we 
will suggest how large cyanides and their corresponding isocyanide isomers may be produced 
in the extended regions surrounding Sgr B2(N). 

The Sgr B2(N) region is known to contain widespread shocks (Chengular & Kanekar 
2003). Such shocks have been presumed to be responsible for the formation and distribution 
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of large aldehyde molecules toward this region (Hollis et al. 2004a, b). Furthermore, it is 
also well known that there is ambient UV flux and cosmic ray flux in these regions. On 
the surfaces and mantles of dust grains, larger species such as CH 3 CN and HC2CN may 
be present. CH 2 CHCN and CH 3 CH 2 CN are formed as a result of H-addition reactions, 
thus depleting HC2CN. As a shock passes, these large molecules are expelled from the grain 
surface into the gas phase. At this point UV and cosmic rays may interact with the cyanide 
species and form isocyanides. As a result, the isocyanides show an extended distribution. In 
hot cores, there is a lower UV and cosmic ray flux due to the large column density of gas 
and dust. Furthermore, thermal processes are the primary mechanism to liberate molecules 
off grain surfaces (Tielens & Hagen 1982; Hasegawa, Herbst, & Lueng 1992). Indeed, this 
is initially why these regions were believed to be the ideal places to form large molecules 
(Ehrenfreund & Charnley 2000). Also, in these compact regions, it is believed warm gas 
phase chemistry in addition to grain surface chemistry can form large molecules which may 
account for the enhanced abundances of CH2CHCN and CH3CH2CN in hot core regions (Liu 
& Snyder 1999, see also CH 3 CH 2 CN column density in Table 2 with respect to Table 1). 
However, thermal processes are not enough to form isocyanide species after the cyanides have 
been expelled from the grain and hence there should be a depletion of isocyanides towards 
hot cores. Our observations support this hypothesis. 

In summary, we have calculated the chemical bonding energies for cyanide molecules 
studied herein (CH 3 CN, CH 2 CHCN, CH 3 CH 2 CN, & HC 4 CN) and found they are much 
greater (>8500 cm" 1 ) than that of the corresponding isocyanide isomers. As a consequence, 
cyanide interstellar molecules are most likely preferentially formed at the expense of iso- 
cyanide isomers. Moreover, in the case of CH 3 CN, it is seen in the LMH hot core and in 
the surrounding medium while CH 3 NC appears to be deficient in the LMH hot core. This 
suggests that large-scale, non-thermal processes such as shocks or enhanced UV flux in the 
surrounding medium account for the conversion of CH 3 CN to CH 3 NC while the LMH hot core 
is dominated by thermal processes which cannot produce a significant amount of CH 3 NC. 
Hudson & Moore (2004) have shown that radiation bombardment of CH 3 CN in ice analog 
experiments can produce CH 3 NC, thus tending to support what is observed. Additionally, 
since CH 3 NC is deficient in the LMH and CH 3 CN is prevalent, this suggests that the relative 
abundance ratio between these two isomers is not meaningful because the molecules are not 
co-spatial. This may well apply to other cyanide and isocyanide isomers. 

J.M.H. gratefully acknowledges research support from H.A. Thronson, Assistant Asso- 
ciate Administrator for Technology, NASA Science Mission Directorate. We also greatly 
thank an anonymous referee whose comments and suggestions provided additional insight to 
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Table 1. GBT Observations - Cyanide and Isocyanide Molecular Line Parameters and Column 

Densities 



Species 


Frequency 


Transition 


S 


E„ 


T c 


(+64 km 


s" 1 ) 


(+73 km 


s" 1 ) 


(+82 km 


s" 1 ) 


N T 


N T (XNC) 














(ATI) 


(AV) 


(ATI) 


(AV) 


(ATI) 


(AV) 


XlO -13 


N T (XCN) 




(MHz) 






(K) 


(K) 


(mK) 


(km/s) 


(mK) 


(km/sj 


(mK) 


(km/s) 


(cm ) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


CH 3 CN a 


18396.7252(7) 


lo-0o 


1.332 


0.88 


10.6 


-275(25) 


5.8(7) 


160(20) 


7.9(9) 


<7 












F = 1 - 1 


























18397.9965(6) 


lo-Oo 


2.220 


0.88 


10.6 


-400(15) 


7.0(3) 


310(19) 


5.9(4) 


<7 




10.8(20) 








F = 2 - 1 


























18399.8924(3) 


lo-0o 


0.444 


0.88 


10.6 


-180(25) 


8.3(9) 


<7 




<7 






0.02(2) 






F = - 1 
























CH 3 NC b 


20105.754(1) 


lo-0o 


4.000 


0.88 


10.2 


-48(2) 


8.4(4) 


11(3) 


6.4(9) 


<7 




0.2(1) 




CH 2 CHCN C 


18512.158(13) 


2 12-lll 


0.375 


3.51 


11.1 


-62(2) 


4.9(1) 


<7 




<7 












F = 2 - 1 


























18513.311(7) 


2 12-lll 


0.701 


3.51 


11.1 


-82(2) 


4.9(1) 


<7 




<7 




37.2(30) 








F = 3 - 2 


























18514.400(27) 


2i2-lll 


0.167 


3.51 


11.1 


-19(3) 


4.1(5) 


<7 




<7 






0.005(2) 






F = 1 - 
























CH 2 CHNC d 


19993.916(1) 


2 12-lll 


1.500 


3.67 


10.5 


8(2) 


1.5(4) 


<7 




<7 




0.2(5) 




CH 3 CH 2 CN e 


18377.721(3) 


2ll-ll0 


1.500 


1.56 


10.0 


-24(2) 


2.8(2) 


20(2) 


2.8(2) 


-20(2) 


2.8(2) 


3.5(4) 






























<0.3 


CH 3 CH 2 NC f 


18803.209(10) 


2i2-lll 


1.500 


1.56 


9.8 


<7 




<7 




<7 




<1.1 




HC 4 CNS 


18638.617(2) 


7-6 


7.000 


2.70 


10.9 


-63(2) 


5.8(2) 


97(1) 


16.0(5) 


-59(2) 


8.0(3) 


6.1(2) 






























<0.03 


HC 4 NC h 


19616.504(2) 


7-6 


7.000 


2.84 


9.7 


<5 




<5 




<5 




<0.2 





Note. — Uncertainties in the frequencies refer to the least significant digit (Taylor & Kuyatt 1994J and are 2 a values (coverage factor 2). Upper limits of the 
line intensities are 3 <r values. All other listed uncertainties arc 1 a. 

a R.cst Frequencies from Boucher et al. 1980; fi a — 3.922(1) from Gadhi et al. 1984 

b Rest Frequencies from Kukolich 1972; Ma=3. 830(60) from Ghosh et al. 1953 

c Rcst Frequencies from Gerry et al. 1979; fi a =3.815(12) , fi b =0. 894(68) from Stolzc, M. & Sutter D. H. 1985 

d Rest Frequencies from Yamada &; Winnewisscr 1975; ^ a =3.470(30) , Mi>=0. 790(100) from Bolton, Owen, & Sheridan 1970 

c Rcst Frequencies, ^ a =3. 850(10) and ^ = 1.230(60) from Lovas 1982 

f Rest Frequencies from Flicgc & Driczlcr 1985; ^a=3. 790(20) and Li b = l. 310(40) from Anderson & Gwinn 1976 
8 Rcst Frequencies from Winnewisscr et al. 1982; /_i a =4.330(30) from Alexander et al. 1976 
h Rcst Frequencies from Botschwina et al. 1998; /j, a =3.590 calculated in this work 
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Table 2. BIMA Array Observations - Molecular Line Parameters and Column Densities 



Species 


Frequency 


Transition 


Sn 2 


E u 


(+64 km 


a" 1 ) 


(+73 km 


s- 1 ) 


(+82 km s _1 ) 


N T 












(Al) 


(AV) 


(A/) 


(AV) 


(A/) 


(AV) 






(MHz) 




(Dcbye 2 ) 


(K) 


(J/bm) 


(km/s) 


(J/bm) 


(km/s) 


(J/bm) 


(km/s) 


(cm- 2 ) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


H 2 CCCC a 


80383.96(2) 


9o9"8oS 


153.8 


19.3 


<0.3 




<0.3 




<0.3 




<6.0xl0 14 


HCOOCH 3 b 


80395.135(126) 


928-909 E 


0.7 


28.8 


<0.3 




<0.3 




<0.3 




<4.7xl0 17 


CH 3 CH 2 CN 


80404.898(12) 


928-827 


126.1 


23.8 


4.6(1) 


16.0(8) 


1.7(3) 


5.1(6) 


0.6(2) 


13.6(4.0) 


3.1 X10 16 


CH3NC 


80405.504(120) 


4 3 -3 3 


29.6 


73.0 


<0.3 




<0.3 




<0.3 






CH 3 NC 


80414.660(120) 


4 2 -3 2 


50.7 


37.4 


<0.3 




<0.3 




<0.3 






CH 3 NC 


80420.062(120) 


4i-3i 


63.3 


16.0 


<0.3 




<0.3 




<0.3 




<4.0xl0 13 


CH 3 NC C 


80421.910(120) 


4 -3o 


67.6 


8.9 


<0.3 




<0.3 




<0.3 





a Rest Frequency and Line Strength (S/J- 2 ) from Killian et al. 1990 
^Rcst Frequency and Line Strength (S^i 2 ) from Bauder 1979 
c Nf calculated using the strongest expected component. 



Table 3. Chemical Bonding 
Energy 



Species 


ZPCBE a 


Relative 






ZPCBE 




(cm" 1 ) 


(cm" 1 ) 


(1) 


(2) 


(3) 


CH3CN 


-29,070,027 





CH3NC 


-29,060,541 


9,486 


CH2CHCN 


-37,408,030 





CH2CHNC 


-37,399,372 


8,658 


CH3CH2CN 


-37,672,368 





CH3CH2NC 


-37,663,671 


8,697 


HC4CN 


-53,818,766 





HC4NC 


-53,807,547 


11,219 



Note. — The Chemical Bonding Energies were 
Obtained at the MP2/aug-cc-pVTZ Level of The- 
ory (Frisch ct al. 2004; Dunning 1989; Kendall, 
Dunning & Harrison 1992; Woon & Dunning, 
1993) 



a Zero-Point Corrected Bonding Energy 
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Fig. 1. — GBT cyanide and isocyanide spectra toward SgrB2(N) at 24.4 kHz channel spacing. 
Transition quantum numbers are shown in each panel. Each spectrum was processed with 
a median filter to remove instrumental slopes in the bandpass. The abscissa is the radial 
velocity with respect to the LSR calculated for the rest frequency of the transition shown 
(see Table 1) at an assumed source velocity of +64 km s _1 . Dashed lines show LSR velocities 
at +64, +73 and +82 km s _1 . In the case of CH 3 CN and CH 2 CHCN the HF structure is 
also seen and labeled for an LSR velocity of +64 km s _1 . 

Fig. 2. — BIMA array spectral passband toward the LMH hot core of Sgr B2(N) containing 
the Jk = 4:r — 3x (K = — 3) transitions of CH 3 NC. In this case, the rest frequency 
corresponds to the CH3CH2CN transition at 80.4049 GHz for a LSR velocity of +64 km s _1 . 
However, there are also emission components at +73 and +82 km s _1 . Finally, the locations 
of the transitions of CH 3 NC, HCOOCH 3 and H 2 CCCC (Table 2) are also labeled assuming 
an LSR velocity of +64 km s _1 . 
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Figure 1. 
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